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Abstract 

The magnitude of the Cabibbo-Kobayashi-Maskawa matrix element 14b bas been measured 
using Y)*^l~v decays recorded on the Z° peak using the OPAL detector at LEP. The D*+ 

D^Tr^ decays were reconstructed both in the particular decay modes D° — » K~7r+ and D° 
K^TT+TT*^ and via an inclusive technique. The product of |T4b| and the decay form factor of the 
go _^ D*+€-i? transition at zero recoil ^(1) was measured to be JP(l)|V;b I = (37.1±1.0±2.0) x IQ-^, 
where the uncertainties are statistical and systematic respectively. By using Heavy Quark Effective 
Theory calculations for T(X), a value of 

jV^bl = (40.7 ± 1.1 ± 2.2 ± 1.6) X 10^3 

was obtained, where the third error is due to theoretical uncertainties in the value of T{V\ The 
branching ratio Br(BO W^^l'v) was also measured to be (5.26 ± 0.20 ± 0.46) %. 
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1 Introduction 

Within the framework of the Standard Model of electroweak interactions, the elements of the Cabibbo- 
Kobayashi-Maskawa mixing matrix are free parameters, constrained only by the requirement that the 
matrix be unitary. The values of the matrix elements can only be determined by experiment. Heavy 
Quark Effective Theory (HQET) provides a means to extract the magnitude of the element Vch from 
particular semileptonic b decays, with relatively small theoretical uncertainties Q. 

In this paper, the value of \ Vch \ is extracted by studying the decay rate for the process B'' — > D*+£~i/ 
as a function of the recoil kinematics of the D*"*" meson [| ||2|-Q. The decay rate is parameterised as 
a function of the variable a;, defined as the scalar product of the four- velocities of the D*"*" and B*^ 
mesons. This is related to the square of the four-momentum transfer from the B*^ to the i'u system, 
g^ by 

_ ml,+ + - 
^~ 0^ — 7^ ' 

and ranges from 1, when the D*^ is produced at rest in the B" rest frame, to about 1.50. Using 
HQET, the differential partial width for this decay is given by 

dT{W^W+l-u) 1 dBr(BO ^ ^,.^2,.^., ,2 ..x 

= = rv.\yj)J- \^)\Vch\ , (2) 

dLO Tgo da; 

where lC{oj) is a known phase space term and ^(w) is the hadronic form factor for this decay |^. 
Although the shape of this form factor is not known, its magnitude at zero recoil, a; = 1, can be 
estimated using HQET. In the heavy quark limit (mb — > oo), ^{uj) coincides with the Isgur-Wise 
function |^ which is normalised to unity at the point of zero recoil. Corrections to J-{'i) have been 
calculated to take into account the effects of finite quark masses and QCD corrections, yielding the 
value and theoretical uncertainty ^(1) = 0.913 ± 0.042 Since the phase space factor lC{w) tends 
to zero as a; — > 1, the decay rate vanishes at u; = 1 and the accuracy of the extrapolation relies on 
achieving a reasonably constant reconstruction efficiency in the region near lo = 1. 

Previous measurements of \Vch\ have been made using B mesons produced on the T{AS) resonance 
and in Z*^ decays 0-^]. These analyses used a linear or constrained quadratic expansion of J-{uj) 
around cj = 1. An improved theoretical analysis, based on dispersive bounds and including higher 
order corrections, has since become available [^0|. This results in a parameterisation for in terms 

of J-{\-) and a single unknown parameter constrained to lie in the range -0.14 < p2 < 1.54^ p2 
corresponding to the slope of !F{uj) at zero recoil. 

The previous OPAL measurement used the decay chain D*^ D^vr^, with the D'^ meson 
being reconstructed in the exclusive decay channels D'^ K^vr"^ and — > K~7r^"7r'^. In this paper, 
a new analysis is described in which only the vr"*" from the D*^ decay is identified, and no attempt 
is made to reconstruct the D*^ decay exclusively. This technique, first employed by DELPHI fs], pT||, 
gives a much larger sample of B*^ — > Y)*^l~v decays than the previous measurement, but also larger 
background, requiring a rather more complex analysis. The measurement of [0] is also updated to use 



the new parameterisation of ^(w) [10|, and improved background models and physics inputs. In both 
cases, the initial number of B" mesons is determined from other measurements of B^ production in TP 
decays. 

The new reconstruction technique is described in Section |2|, the determination of uj for each event 
in Section the fit to extract J-'(l)|V^b| ^'iid in Section |^ and the systematic errors in Section H. 
The updated exclusive measurement is discussed in Section ^ and the measurements are combined 
and conclusions drawn in Section 0. 



^Charge conjugate reactions are always implied, and the symbol £ refers to either an electron or muon. 
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2 Inclusive reconstruction of D*+£-i> events 



The OPAL detector is well described elsewhere The data sample used in this analysis consists of 
about 4 million hadronic decays collected during the period 1991-1995, at centre-of-mass energies 
in the vicinity of the resonance. Corresponding simulated event samples were generated using 
JETSET 7.4 ill as described in |ll. 



Hadronic Z'^ decays were selected using standard criteria [|l^. To ensure the event was well 
contained within the acceptance of the detector, the thrust axis direction^ was required to satisfy 
\ cos 9t\ < 0.9. Charged tracks and electromagnetic calorimeter clusters with no associated tracks 



were then combined into jets using a cone algorithm |15], with a cone half angle of 0.65 rad and a 
minimum jet energy of 5 GeV. The transverse momentum pt of each track was defined relative to the 
axis of the jet containing it, where the jet axis was calculated including the momentum of the track. 
A total of 3 117 544 events passed the event selection. 

The reconstruction of B'^ — >■ D*~^£~i' events was then performed by combining high p and pt lepton 
(electron or muon) candidates with oppositely charged pions from the D*'^ — > D^vr^ decay. Electrons 
were identified and photon conversions rejected using neural network algorithms |jl^ , and muons were 
identified as in [16|. Both electrons and muons were required to have momenta p > 2 GeV, transverse 



momenta with respect to the jet axis pt > 0.7 GeV, and to lie in the polar angle region | cos6\ < 0.9. 

The event sample was further enhanced in semileptonic b decays by requiring a separated secondary 
vertex with decay length significance L/ai > 2 in any jet of the event. The vertex reconstruction 



algorithm and decay length significance calculation are described fully in |14|. Together with the 
lepton selection, these requirements result in a sample which is about 90 % pure in bb events. 

An attempt was made to estimate the direction in each jet containing a lepton candidate. Each 
track (apart from the lepton) and calorimeter cluster in the jet was assigned a weight corresponding 
to the estimated probability that it came from the D'' decay. The track weight was calculated from 
an artificial neural network, trained to separate tracks from b decays and fragmentation tracks in 



b jets 1 14]. The network inputs are the track momentum, transverse momentum with respect to the 
jet axis, and impact parameter significances with respect to the reconstructed primary and secondary 
vertices (if existing). The cluster weights were calculated using their energies and angles with respect 
to the jet axis alone, the energies first being corrected by subtracting the energy of any charged tracks 



associated to the cluster [|17|. 



Beginning with the track or cluster with the largest weight, tracks and clusters were then grouped 
together until the invariant mass of the group (assigning tracks the pion mass and clusters zero mass) 
exceeded the charm hadron mass, taken to be 1.8 GeV. If the final invariant mass exceeded 2.3 GeV, 
the jet was rejected, since Monte Carlo studies showed such high mass candidates were primarily 
background. For surviving jets, the momentum p^o of the group was used as an estimate of the D*^ 
direction, giving RMS angular resolutions of about 45 mrad in (j) and 0. The D*^ energy was calculated 
as E-QO = 

The selection of pions from D*^ decays relies on the small mass difference of only 145 MeV [^] 
between the D*"*" and D^, which means the pions have very little transverse momentum with respect 
to the D'^ direction. Each track in the jet (other than the lepton) was considered as a slow pion 
candidate, provided it satisfied 0.5 < p < 2.5 GeV and had a transverse momentum with respect to 
the D'^ direction of less than 0.3 GeV. If the pion under consideration was included in the reconstructed 
DO, it was removed and the momentum and energy recalculated. The final selection was made 
using the reconstructed mass difference Am between D*"*" and mesons, calculated as 



Am = \J E"^, - Ipd* P - moo , 

where the D*^ energy is given by E-q* = E-^o + Et^ and momentum by pD* = Pd" + Ptt- 

right handed coordinate system is used, with positive z along the electron beam direction and x pointing to the 
centre of the LEP ring. The polar and azimuthal angles are denoted by 9 and cj>. 
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A new secondary vertex was then iteratively reconstructed around an initial seed vertex formed 
by the intersection of the lepton and slow pion tracks. Every other track in the jet was added in turn 
to the seed vertex, and the vertex refitted. The track resulting in the lowest vertex fit ^ was retained 
in the seed vertex for the next iteration. The procedure was repeated until no more tracks could be 
added without reducing the vertex fit probability to less than 1%. The decay length L' between 
the primary vertex and this secondary vertex, and the associated error 0"^', were calculated as in 
The pion candidate was accepted if the decay length satisfied —0.1cm < L' < 2 cm and the decay 
length significance satisfied L' joi^i > —2. 

The resulting distributions of Am for opposite and same sign lepton-pion combinations are shown 
in Figure [l|( a) and (b). The predictions of the Monte Carlo simulation are also shown, broken down 
into contributions from signal B'' — > D*~^i~i' events, 'resonant' background containing real leptons and 
slow pions from D*"*" decay, and combinatorial background, made up of events with fake slow pions, 
fake leptons or both. 

In Monte Carlo simulation, about 45% of opposite sign events with Am < 0.17 GeV are signal 
BO ^ D*+£-i? events, 14 % are resonant background and 41 % are combinatorial background. The 
resonant background is made up mainly of B~ — > D*~^iT~i~iy, — > D*~^iT^i~i' and Bg D*~^K^i~i' 
decays. These are expected to be dominated by b semileptonic decays involving orbitally excited 
charm mesons (generically referred to as D**), e.g. B~ — > D**^£~i> followed by D**° — > D*"'"7r~. 
These decays will be denoted collectively by B — > D*^hi~v'. Small contributions are also expected 
from b D*"'"tPX decays with the r decaying leptonically, and b — »■ D*+D~X with the D~ decaying 
semileptonically (each about 1 % of opposite sign events). For same sign events with Am < 0.17 GeV, 
there is a small resonant contribution of about 6 % from events with a real D*"*" — > D'^vr^ where the 
D*^ decays semileptonically, and the rest is combinatorial background. 

The most important background, from B D*~^h.i~i' decays, comes from both charged B+ and 
neutral B'^ and Bg decays, whereas the signal comes only from B'' decays. Therefore the charge Qvtx of 
the reconstructed secondary vertex containing the lepton and slow pion and its estimated error ctq^^^ 
were calculated, using 

i 

^Wi{l - Wi)qf , 

i 

where Wi is the weight for track i of charge qi to come from the secondary vertex, and the sums are 
taken over all tracks in the jet []l9| , p0| . The weights were calculated in a similar way to those used for the 
D*^ direction reconstruction, using a neural network with the track momentum, transverse momentum 
and impact parameter significances with respect to the reconstructed primary and secondary vertices 
as inputs. The weights for the lepton and slow pion candidate tracks were set to one. The vertex 
charge distributions for opposite sign events with Am < 0.17 GeV are shown in Figure |l|(c) and (d). 

The reconstructed vertex charge and error were used to divide the data into different classes 
enhanced or depleted in B+ decays, thus reducing the effect of this background and increasing the 
statistical sensitivity. Three classes c were used — class 1 where the charge is measured poorly (ctq^^,^ > 
0.9), class 2 where the charge is measured well and is compatible with a neutral vertex {ctq^^^ < 0.9, 
IQvtxl < 0.5) and class 3 where the charge is measured well and is compatible with a charged vertex 
(f^Qvtx < 0.9, IQvtxl > 0.5). 

3 Reconstruction of u: 

The recoil variable lo was estimated in each event using the reconstructed four-momentum transfer to 
the iu system: 

g2 = (EB0-SD.)'-(PB0-PD*)' 



Q 



vtx 



cr. 



Qvtx 
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vertex charge vertex charge Q, 



Figure 1: Reconstructed Am distributions for selected (a) opposite sign and (b) same sign lepton-pion 
combinations; reconstructed vertex charge Qvtx distributions for opposite sign events with Am < 
0.17 GeV and (c) ctq^.j^ > 0.9 and (d) ctq^^.^ < 0.9. In each case the data are shown by the points 
with error bars, and the Monte Carlo simulation contributions from signal B'' Y)*'^i^u decays, other 
resonant D*"*" decays and combinatorial background are shown by the open, single and cross hatched 
histograms respectively. The vertex charge classes 1 (poorly reconstructed), 2 (well reconstructed 
neutral vertex) and 3 (well reconstructed charged vertex) are indicated. 
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together with equation |T|. The B'' and D*"*" energies £"30 and -Ed* were estimated directly, whilst the 
momentum vectors pgo and po* were estimated using the energies together with the reconstructed 
polar and azimuthal angles, as described in more detail below. 

Since the slow pion has very little momentum in the rest frame of the decaying D*"^, the momentum 
(and hence energy) of the D*"*" was estimated by scaling the reconstructed slow pion momentum by 
m-£,*+ /m-,^, as for the D*^ — > K~7r''~7r'' channel in [^. A small correction (never exceeding 12%) was 
applied, as a function of cos 9*, the angle of the slow pion in the rest frame of the D^. This procedure 
gave a fractional D*"*" energy resolution of 15%. The polar and azimuthal angles of the D*"*" were 
reconstructed using weighted averages of the slow pion and D*^ directions, giving resolutions of about 
22 mrad on both (p and 9. 

The energy of the B*^ was estimated using a technique similar to that described in ||2^, exploiting the 
overall energy and momentum conservation in the event to calculate the energy of the unreconstructed 
neutrino. First, the energy -Ebjet of the jet containing the B" was inferred from the measured particles 
in the rest of the event, by treating the event as a two-body decay of a into a b jet (whose mass was 
approximated by the B" mass) and another object making up the rest of the event. Then, the total 
fragmentation energy .Bfrag in the b jet was estimated from the measured visible energy in the b jet 
and the identified B*^ decay products: -Efrag = E^is — — Ed* . Finally, the B*^ energy was calculated 
as E-QO = -Ebjet — -Efrag) giving an RMS resolution of 4.4 GeV. 

The b direction was estimated using a combination of two techniques. In the first, the momentum 
vector of the B*^ was reconstructed from its decay products: 

Pbo = Pd* +Pe + Pu, 

the neutrino energy being estimated from the reverse of the event visible momentum vector: — 
— Pvis- The visible momentum was calculated using all the reconstructed tracks and clusters in the 
event, with a correction for charged particles measured both in the tracking detectors and calorimeters 
[^. This direction estimate is strongly degraded if a second neutrino is present {e.g. from another 
semileptonic b decay in the opposite hemisphere of the event), and its error was parameterised as a 
function of the visible energy in the opposite hemisphere. The resulting resolution is typically between 
40 and 100 mrad for both (p and 9. 

The second method of estimating the B*^ direction used the vertex flight direction — i.e. the di- 
rection of the vector between the primary vertex and the secondary vertex reconstructed around the 
lepton and slow pion as described in Section |^. The accuracy of this estimate is strongly dependent on 
the B*^ decay length, and was used only if the decay length significance L' /a^i exceeded 3. After this 
cut, the angular resolution varies between about 15 and 100 mrad, and is worse for 9 in the 1991 and 
1992 data where accurate z information from the silicon microvertex detector was not available. The 
go B*+i-u candidate was rejected if the two reconstruction methods gave 9 01 (f) angles disagreeing 
by more than three standard deviations, which happened in 7 % of Monte Carlo signal events. Finally, 
the two B*^ direction estimates were combined according to their estimated uncertainties, giving aver- 
age resolutions of 35 mrad on (p and 43 mrad on 9, including events where only the first method was 
used. 

The estimate of derived from the B*^ and D*^ energies and angles was improved by applying 
the constraint that the mass /i of the neutrino produced in the B*^ T)*^l^9 decay should be zero. 
The neutrino mass is given from the reconstructed quantities by: 

/i^ = {E^o - Ed* - Eif - (pbo - PD- - Pif . 

The constraint was implemented by calculating using a kinematic fit, incorporating the measured 
values and estimated uncertainties of the B'' and D*"*" energies and angles (the B'^ angular uncer- 
tainties varying event by event). This procedure improved the average resolution from 2.78 GeV^ 
to 2.57 GeV^. In Monte Carlo simulation, 11% of signal events were reconstructed with < 0, 
corresponding to an unphysical value of uj larger than Wmax ~ 1-5, and were rejected. 
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The resulting distributions of reconstructed u! for various ranges of true lo (denoted to') in simulated 
BO B*+i-u decays are shown in Figure |2|(a-e). The average RMS uj resolution is about 0.12, but 
there are significant non-Gaussian tails. The resolution was parameterised (separately for the 1991-2 
and 1993-5 data) as a continuous function R{ijj,uj') giving the expected distribution of reconstructed 
LO for each true value oj' . The resolution function was implemented as the sum of two asymmetric 
Gaussians (i.e. with different widths either side of the peak) whose parameters were allowed to vary 
as a function of oo' . The convolution of this resolution function with the Monte Carlo lo' distribution 
is also shown in Figure ||(a-e) , demonstrating that the resolution function models the lo distributions 
well. 

The efficiency to reconstruct B*^ — > D*^^~i> decays is shown in Figure |2|(f), together with a 
second order polynomial parameterisation. The efficiency varies with lo' , but is reasonably flat in the 
critical region near oj' = 1 where the extrapolation to measure ^(l)|V^b| is carried out. 

4 Fit and results 

The values of ^(l)|V^b| ^-nd were extracted using an extended maximum likelihood fit to the 
reconstructed mass difference Am, recoil lo and vertex charge class c of each event. Both opposite and 
same sign events with Am < 0.3 GeV and uo < t^max were used in the fit, the high Am and same sign 
events serving to constrain the combinatorial background normalisation and shapes in the opposite 
sign low Am region populated by the B*^ decays. Using the Am value from each event in 

the fit, rather than just dividing the data into low Am 'signal' and high Am 'sideband' mass regions, 
increases the statistical sensitivity as the signal purity varies considerably within the low Am region. 
The logarithm of the overall likelihood was given by 

ln/: = ^ln/:^ + ^ln4-A^'^-iV'' (3) 
i=i j=i 

where the sums of individual event log-likelihoods ln.Cf and ln£^ are taken over all the observed 
opposite sign and M'^ same sign events in the data sample, and N"" and A^'' are the corresponding 
expected numbers of events. 

The likelihood for each opposite sign event was given in terms of different types or sources of event 

by 

4 

£^(Am„^„Q) = Y.N^fl,Ms,cA^m,)Ps{Lo,) (4) 

where A'^^ is the number of expected events for source s, /"^ is the fraction of events in source s 
appearing in vertex charge class c, Ms^c(^'tt-) is the mass difference distribution for source s in class 
c and -Ps(w) the recoil distribution for source s. For each source, the vertex charge fractions /^^ sum 
to one and the mass difference Mg^d^'tTT'i) and recoil Ps{tjj) distributions are normalised to one. The 
total number of expected events is given by the sum of the individual contributions: A^" = J2t=i ^s- 
There are four opposite sign sources: (1) signal B" — > events, (2) B D*^hi~v' events 

where the D*+ is produced via an intermediate resonance (D**), (3) other opposite sign background 
involving a genuine lepton and a slow pion from D*"*" decay and (4) combinatorial background. The 
sum of sources 2 and 3 are shown as 'resonant background' in Figure |l|. A similar expression to 
equation ^ was used for the event likelihood for same sign events. In this case, only sources 3 
and 4 contribute. 

The mass difference distributions Mc^s(Am) for sources 1-3 were represented by analytic functions, 
whose parameters were determined using large numbers of simulated events, as were the recoil distri- 
butions Ps{uj) for sources 2 and 3. The fractions in each vertex charge class for sources 1-3 were also 
taken from simulation. For the signal (source 1), the product of the expected number of events 
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Figure 2: Reconstructed distributions for various ranges of true to (denoted lo') in Monte Carlo 
B'' — > D*+^~P events, together with the prediction from the resolution function. Events reconstructed 
in the unphysical region with uj > o^max are rejected. The reconstruction efficiency as a function of 
true Lo is also shown. 
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Quantity 


Assumed value 


Reference 


Rh 

/bo 
Tbo 

Br(D*+ ^ D07r+) 
Br(b ^ D*+heu) 
Br(b ^ D*+r"PX) 
Br(BO ^D*+dW") 
Br(b ^ D*+X) 
Br(c ^ D*+X) 
Br(DO ^ K-7r+) 
Br(DO ^ K~7r+7r0) 


(21.70 ± 0.09) % 
(39.7i^|) % 
1.56 ± 0.04 ps 
(68.3 ± 1.4)% 

(0.76 ±0.16)% 

(0.65 ±0.13)% 
(4.2 ± 1.5) % 
(17.3 ± 2.0) % 
(22.2 ± 2.0) % 

(3.85 ± 0.09) % 
(13.9 ± 0.9) % 


I 
[ 
[ 
[ 

I 
[ 
1 
1 

I 
I 


18 
18 
18 
18 
22 
18 
18 
24 
24 
18 
18 


1 
1 
1 
1 

, see text 
, see text 

1 
1 
1 

1 
1 



Table 1: Input quantities used in the fits for jr(l)|yc.b| and p^. The values marked 'see text' are 
derived using methods explained in Section ^ 



and recoil distribution -Pi(w) was given by convolving the differential partial decay width (equation ^) 
with the signal resolution function and reconstruction efficiency: 

J^^{uj')\Vcb\'^!C{L0')e{L0')R{L0,Uj')di0' 

where Nz is the number of hadronic decays passing the event selection, i?b = r^^f^/Fhad is the 
fraction of hadronic decays to bb, /go the fraction of b quarks hadronising to a B*^ and rgo the 
B^ lifetime. The factor of four accounts for the two b hadrons produced per — > bb event and the 



two identified lepton species (electrons and muons). The form factor J-{uj') is given in |10| in terms of 
the normalisation .F(l) and slope parameter p^. The efficiency function e(u;') and resolution function 
R{uj,uj') were described in Section ^, and the known phase space factor }C{ui') is given in The 
assumed values of the numerical quantities are given in Table |^. 

Since the data are divided into different vertex charge classes enhanced and depleted in B"*" decays, 
the fit gives some information on the amount of B ^ T)*^h£~i' background. The predicted level of 
this background in the fit was therefore allowed to vary under a Gaussian constraint corresponding 
to the branching ratio of Br(6 — > T)*~^hiv') = (0.76 ± 0.16)%. The latter has been calculated from 



the measured branching ratio Br(b — > D*"'"7r~£i/X) = (0.473 ± 0.095)% [22|, assuming isospin and 
SU(3) flavour symmetry to obtain the corresponding b — s- D*~^'ir^ii' and b — s- T)*~^K^ii' branching 
ratios. A scaling factor of 0.75 ± 0.25 was included for the last branching ratio to account for possible 
SU(3) violation effects reducing the branching ratio D**"*" — > D*"'"K'^ compared to the expectation of 
|Br(D**^ — > D*"'"7r'^). The P2{u>) distribution for these events was taken from simulation, using the 
calculation of Leibovich et al. l23] to predict their recoil spectrum. 



The numbers N!^'^ and ^3(0;) distributions for the small background contributions covered by 
source 3 (both opposite and same sign) were taken from Monte Carlo simulation, with branching 
ratios adjusted to the values given in Table |l|. The branching ratio for b — > D*"^r~P was derived 
using the inclusive branching ratio Br(b — > rX) = (2.6 ± 0.4) % ||T^, assuming a D*"*" is produced in 
a fraction 0.25 ± 0.03 of the time, as seen for the corresponding decays b — > D*"'"£X and b — > i'K 
{£ = e, /i) IM]. The rate of b ^ D*^" D~X was assumed to be dominated by the two body decay 
B*^ — > D*^Ds* . The rate of real D*^ decays in the same sign background depends on the production 



fractions of D*"*" in bb and cc events, which were taken from |24|. 

The parameters of the analytic functions describing the combinatorial background {N^, N^, f^^, 
f^^, M^^d^iT^) and P4{uj)) were fitted entirely from the data, with only the choice of functional forms 
motivated by simulation. The shapes of the mass and recoil functions (including a small correlation 
between Am and uj) are constrained by the same sign sample (which is almost entirely combinatorial 
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background), and are the same for each vertex charge class c. The opposite sign high Am region 
serves to normalise the number of combinatorial background events in the low Am region for each 
vertex charge class. 

The values of ^(l)|T4b| and were extracted by maximising the total likelihood given by equa- 
tion p. The values of ^(l)|V^b| and p"^ were allowed to vary, together with the level of B ^ D*^li£~i> 
background and 13 auxiliary parameters describing the combinatorial background distributions. A 
result of 

H'^Wcbl = (37.5 ± 1.2) X 10-3, 
p2 = 1.12 ±0.14 

was obtained, where the errors are only statistical. The correlation between .F(l)|V^b| and p'^ is 0.77. 
The distributions of reconstructed to for opposite and same sign events with Am < 0.17 GeV, together 
with the fit results, are shown in Figure ^. The fit describes the data well, both in this region and the 
high Am region dominated by combinatorial background, for all three of the vertex charge classes. 

By integrating the differential partial decay width (equation |2|) over all values of uj, the branching 
ratio Br(B° ^ D*+t u) was also determined to be 

Br(BO ^ D*+r z^) = (5.92 ± 0.27) % , 

where again the error is statistical only. This result is consistent with the world average once systematic 
errors are included. 

Many previous results have been obtained using a constrained quadratic expansion for the form 
factor: = -?^(1) [1 — a?{<jJ — 1) + hiuj — 1)^] , where a is a slope parameter to be determined by the 

fit, and h is constrained to 6 = 0.66a^ — 0.11 To allow comparison with such measurements, the fit 
was also performed with this parameterisation of J'ioj), giving the results .7-'(l)| V^bl = (36.9ibl.2) x 10"^ 
and = 0.88 it 0.14, the correlation between ^(l)|14b| and being 0.79. The difference in the two 
curvature parameters p^ and is in good agreement with the expectation of p^ — a?' ^ 0.21 [10|. 



5 Systematic Errors 

Systematic errors arise from the uncertainties in the fit input parameters given in Table the Monte 
Carlo modelling of the signal u; resolution, the recoil spectrum of b ^ T)**ii/ decays and selection 
efficiencies, and possible biases in the fitting method. The resulting systematic errors on the values of 
•^(l)l^bl) p'^ and Br(B'^ D*~^£~r') are summarized in Table |2| and described in more detail below. 

Input quantities: The various numerical fit inputs were each varied according to the errors given in 
Table || and the fit repeated to assess the resulting uncertainties. 

b — > D**£i5 decays: The calculation of Leibovich et al. was used to simulate the recoil spectrum 
of B — > D*'^h^-i> decays, assumed to be produced via the semileptonic decay B'' — > T)**iu. 
Here D** represents a P-wave orbitally excited charm meson. The calculation predicts the recoil 
spectra and relative rates of semileptonic decays involving both the narrow Di and D2 states 
and wide Dq and DJ states. All these decays are suppressed close to w = 1 by an extra factor 
of {uP' — 1) when compared with the signal B'' — > 'D*~^l~i> decays. This reduces the uncertainty 
due to the rate of B — > D*"'"h£-i^ decays in the extrapolation of the signal recoil spectrum to 
Lo = \. Non-resonant B D*^h£-z/ decays are not included in the model, but are not expected 
to contribute close to tj = 1. 

The differential decay rates in ||2^ are given in terms of five possible expansion schemes and 
several unknown parameters: a kinetic energy term ryke and the slopes of the Isgur-Wise functions 
for the narrow and wide D** states fi and ^1. These parameters were varied within the allowed 
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Figure 3: Distributions of reconstructed uj for (a) opposite sign and (b) same sign events with Am < 
0.17 GeV. The data are shown by the points with error bars and the expectation from the fit result by 
the histograms. The contributions from signal B° — > D*'^£~i/, resonant and combinatorial backgrounds 
are indicated. 
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Error Source 


A{m\Vc^\)/im\Vcb\) (%) 


Ap2 


ABr/Br (%) 


Rh value 


0.2 


- 


0.4 


B'^ lifetime 


1.3 


- 


2.6 


/go value 


+2.8 
-2.3 


- 


+5.5 
-4.5 


Br(D*+ ^ D07r+) 


1.0 


- 


2.1 


Br(b ^ D*+T-u) 


0.1 


0.00 


0.4 


Br(b ^ D*+Di*^") 


0.3 


0.01 


1.3 


Br(b ^ D*+X) 


0.0 


0.00 


0.1 


Br(c ^ D*+X) 


0.0 


0.00 


0.0 


b ^ D**^z^X rate 


0.6 


0.05 


2.9 


b ^ D**£i?X model 


4.1 


0.19 


3.5 


b fragmentation 


1.1 


0.12 


3.2 


D'^ decay multiplicity 


1.2 


0.08 


2.2 


a; reconstruction 


2.2 


0.12 


2.1 


Lepton identification efficiency 


1.2 


- 


2.3 


Vertex tag efficiency 


1.7 




3.4 


Track reconstruction efficiency 


0.1 


0.10 


4.3 


Tracking resolution 


2.4 


0.04 


4.4 


Fitting method 


1.9 


0.05 


0.6 


Total 


6.8 


0.29 


11.5 



Table 2: Summary of systematic errors on the measured values of ^(1) | V^b I > ^-nd Br(B^ — > D*~^i~r') 
for the inclusive analysis. The fractional errors on ^(l)|14b| and Br(B'' — > are given, whereas 

the errors on are absolute. 



ranges —0.75 < r/tc < 0.75 GeV, — 2 < ti < —1 and — 2 < Ci < subject to the constraint that 
the ratio R = r{B D^eiy)/{T{B Diiv) lie within the measured range R = 0.37 ± 0.16 



1 22, 26]. This excludes the expansion schemes A^o and of [23| and constrains the allowable 
values of in the others. The fraction of B ^ D**£~^' decays involving the narrow Di and 
D2 states, which is not precisely predicted, was varied in the range 0.22 it 0.06, obtained by 
comparing the measured rates for B*^ semileptonic decays involving D"*", D*"*", Di and D2 with 
the inclusive semileptonic decay rate p!8|, p^, 26 1. The systematic errors were determined as half 



the difference between the two parameter sets giving the most extreme variations in ^(l)|V^b| 
and p^, and the central values were adjusted to half way between these two extremes. The 
values of both J^(l)|V^b| and are most sensitive to variations in r/ke, which is constrained by 
the measured value of R. 

b fragmentation: The effect of uncertainties in the average b hadron energy {xe) = -E-bZ-E'beam 
was assessed in Monte Carlo simulation by reweighting the events so as to vary {xe) in the 
range 0.702 it 0.008 |2^, and repeating the fit. The largest of the variations observed using the 
fragmentation functions of Peterson, Collins and Spiller, Kartvelishvili and the Lund group [^] 
were taken as systematic errors. 

D° decay multiplicities: The signal reconstruction efficiency and vertex charge distributions are 
sensitive to the B'' decay multiplicity, which depends only on the D'' decay for the B'' decay 
channels of interest. The systematic error was assessed by varying separately the D'' charged 
and vr'' decay multiplicities in Monte Carlo simulation according to the measurements of Mark 



III The branching ratio 0° -> K°,K° was also varied according to its uncertainty [|l^ . 

The resulting uncertainties on ^(l)|ycb| and p^ from each variation were added in quadrature 
to determine the total systematic errors. 
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to reconstruction: The modelling of the uj resolution depends on the description of the D*+ and 
B'' energy and angular distributions in the simulation. The reconstructed D*"*" and energy 
distributions in data and simulation were compared, and the means were found to differ by 0.04 
and 0.13 GeV respectively. The opposite hemisphere missing energy was found to agree within 
5%. The corresponding systematics were assessed by shifting or scaling the data distributions 
and repeating the uj reconstruction and fit. The modelling of the angular resolution was checked 
by studying the agreement of the two angular estimators — the slow pion and D" directions for 
the the D*"*", and the missing energy vector and vertex flight directions for the B'^. The angular 
resolutions were found to be up to 5 % worse in data, and the systematic error was assessed 
by degrading the simulated resolution appropriately. Finally, the fraction of events with lv 
reconstructed in the physical region (to < Wmax) was found to be 3.5 % smaller in the data, in 
both opposite and same sign charge samples. The reconstruction efficiency was corrected for 
this effect, and an additional systematic error of half the correction (1.7%) assumed. The final 
systematic errors due to to resolution modelling are dominated by the B'^ 6 resolution. 

Lepton identification efficiency: The electron identification efficiency has been studied using con- 
trol samples of pure electrons from e^e~ — > e^e~ events and photon conversions, and found to be 



modelled to a a precision of 4.1 % |14]. The muon identification efficiency has been studied using 



muon pairs producted in two-photon collisions and n'^n events, giving an uncertainty of 



2.1% l30|l. 



Vertex tag efficiency: The fraction of hemispheres with identified leptons which also had a selected 
secondary vertex was found to be about 4 % less in data than in simulation. The overall fraction 
of vertex tagged hemispheres was also found to be about 4 % lower in data. These discrepancies 
were translated into systematic errors on the efficiency to tag a semileptonic b decay with a 
secondary vertex in either the same or the opposite hemisphere, in each case attributing the 
whole discrepancy to a mismodelling of b hadron decays. The resulting errors on the same and 
opposite hemisphere tagging efficiencies were taken to be fully correlated. 

Track reconstruction efficiency: The overall track reconstruction efficiency is known to be mod- 



elled to a precision of 1 % |14], and a similar uncertainty was found to be appropriate for the 
particular class of slow pion tracks from D*+ decays. The systematic error was assessed by 
randomly removing 1 % of tracks in the simulation and repeating the fit. 

Tracking resolution: Uncertainties in the tracking detector resolution affect the efficiency, uj recon- 
struction and vertex charge distributions. The associated error was assessed in the simulation 
by applying a global 10 % degradation to all tracks, independently in the r-cj) and r-z planes, as 
in [Hi. 



Fit method: The entire fitting procedure was tested on a fully simulated Monte Carlo sample seven 
times bigger than the data, with true values of J-'(l)|V^b| = 32.5 x 10~^ and = 1.3. The fit 
gave the results J^(l)|V;b| = (31.8 ± 0.5) x lO^^ and = 1.25 ± 0.07, consistent with the true 
values. For each variable, the larger of the deviation of the result from the true value and the 
statistical error were taken as a systematic errors due to possible biases in the fit. Additionally, 
the large Monte Carlo sample was reweighted to change the values of .7^(l)|V^b| and p^, and 
the fit correctly recovered the modified values. To verify the correctness of the statistical errors 
returned by the fit, it was performed on many separate subsamples, and the distribution of fitted 
results studied. Further checks on the data included performing the analysis separately for B" 
decays involving electrons and muons, dividing the sample according to the year of data taking, 
and varying the lepton transverse momentum cut. In all cases, consistent results were obtained. 

Including all systematic uncertaintities, the final result of the inclusive analysis is 

T{l)\Vch\ = (37.5 ± 1.2 ± 2.5) X 10-^ 
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p2 = 1.12 ±0.14 ±0.29, 
where the first error is statistical and the second systematic in each case. 

6 Measurement using exclusively reconstructed D*+ decays 

In this analysis, the D*^ from the B*^ —^ D*^£^ij, D*^ — > D'^vr^ decay is reconstructed exclusively 
in the decay modes D" —> K^vr"*" ('3-prong') and K~7r"'"7r'^ ('satellite' — where the tt" is not 

reconstructed). The event selection, reconstruction and determination of uj are exactly the same as 
described in The determination of the signal and background fractions, and the fit to extract 
•^(l)l^bl and are performed as in but have been updated using the improved form factor 
parameterisation |jlO|, the updated input parameters given in Table |^ and the b — >■ D**ii' decay 
model [23 1 discussed in Section ^. 



The fit is performed on 814 3-prong and 1396 satellite candidates, of which 505 ± 44 and 754 ± 72 
are attributed to B*^ T>*~^£~i' signal decays. The result of the fit is 



T{l)\Vch\ = (36.8 ± 1.6 ± 2.0) X 10~^ 
p2 = 1.31 ±0.21 ±0.16, 

where again the first errors are statistical and the second systematic. The statistical correlation 
between .?^(l)|ycb[ and p'^ is 0.95. The distribution of reconstructed w for selected candidates (both 
3-prong and satellite) is shown in Figure ^. The branching ratio Br(B'^ T)*^£~i>) has also been 
determined to be 



Br(B° = (5.11 ± 0.19 ± 0.49) 



'0. 



The systematic errors arise from uncertainties in the background levels in the selected samples, 
as well as uncertainties in the Monte Carlo simulations. They have been evaluated following the 
procedures described in |^, and are summarised in Table ^. The selection efficiency error includes 
contributions from lepton identification efficiency, b fragmentation and detector resolution uncer- 



tanties, described in detail in |31]. The largest change with respect to the previous result is due to 
the improved b — > D**£v' modelling, with the suppression of this background at values of uj close to 
one. This reduces the statistical error, the systematic uncertainty due to the rate of such decays, and 
shifts the central value of .F(l)|ycb| upwards as compared to p. 



7 Conclusions 

The CKM matrix element |V^b| has been measured by studying the rate of the semileptonic decay 
u as a function of the recoil kinematics of both inclusively and exclusively reconstructed 
D*+ mesons. The two results are combined, taking into account the statistical correlation of 18 % and 
correlated systematic errors from physics inputs and detector resolution. The results are: 

J^{l)\Vcb\ = (37.1 ± 1.0 ±2.0) X 10~^ 
p2 = 1.21 ±0.12 ±0.20, 
Br(BO ^ D*+i-u) = (5.26 ± 0.20 ± 0.46) % , 

where the first result is statistical and the second systematic in each case. The statistical and sys- 
tematic correlations between J^(l)|14b| and p'^ are 0.90 and 0.54 respectively. These results supersede 
our previous publication [^. They are consistent with other determinations of ^(l)|V^b| at LEP ||^^ 
and the T(4S') resonance [§. The branching ratio is consistent with the world average result of 
(4.60 ± 0.27) % [|l8|. The result for .F(l)|ycb| is the most precise to date from any single experiment. 
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Figure 4: Distribution of reconstructed lo for B*^ — > 'D*~^l~v candidates in the exclusive analysis. The 
data are shown by the points with statistical error bars, and the fit result by the histogram. The 
contributions from signal B'' resonant and combinatorial backgrounds are indicated. 
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Error Source 


A{T{l)\V,b\)/{T{l)\Vcb\) (%) 


Ap2 


ABr/Br (%) 


Rh value 


0.2 




0.4 


B" lifetime 


1.2 


_ 


2.4 


/go value 


2.5 


- 


4.9 


Br(D*+ ^ D07r+) 


1.0 


- 


2.0 


Br(DO ^ K-7r+) 


0.4 


0.011 


0.2 


Br(D° ^ K-7r+7r°) 


2.6 


0.025 


3.8 


Br(b ^ D*+r-p) 


0.1 


- 


0.2 


Br(b^D*+Di*^") 


0.2 




0.4 


b ^ D**£z^X rate 


1.1 


0.078 


2.0 


b ^ D**^z^X model 


1.0 


0.128 


5.0 


Combinatorial background 


1.2 


0.012 


1.7 


Fake leptons 


0.2 




0.4 


Fake D" 


0.5 


0.005 


0.7 


w reconstruction 


1.4 


0.035 




Selection efficiency 


2.9 


0.005 


3.1 


Total 


5.5 


0.157 


9.5 



Table 3: Summary of systematic errors on the measured values of .7-"(l)|V^b|) ^-nd Br(B'^ — s- D*^"£~z/) 
for the exclusive analysis. The fractional errors on ^(l)|V^b| ^-^d Br(B'^ — > D*~^£~i') are given, whereas 
the errors on are absolute. 

Using the theoretical estimate ^{1) = 0.913 ± 0.042 j^, the value of \Vch\ is determined to be 

= (40.7 ± 1.1 ± 2.2 ± 1.6) X 10"^ 
where the uncertainties are statistical, systematic and theoretical respectively. 
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